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Ph Abstract 

^ The dynamic initiation of sliding at planar interfaces between deformable and rigid solids 

O 

c/5 is studied with particular focus on the speed of the slip front. Recent experimental results 

O 

' ^ showed a close relation between this speed and the local ratio of shear to normal stress 

>-. 

^ measured before slip occurs (static stress ratio). Using a two-dimensional finite element 

'""' model, we demonstrate, however, that fronts propagating in different directions do not 

^ have the same dynamics under similar stress conditions. A lack of correlation is also 

^ observed between accelerating and decelerating slip fronts. These effects cannot be entirely 

O 

\0 associated with static local stresses but call for a dynamic description. Considering a 

Q dynamic stress ratio (measured in front of the slip tip) instead of a static one reduces the 

,-H above-mentioned inconsistencies. However, the effects of the direction and acceleration are 

. ^ still present. To overcome this we propose an energetic criterion that uniquely associates, 

^ independently on the direction of propagation and its acceleration, the slip front velocity 
d 

with the relative rise of the energy density at the slip tip. 



Introduction 

Many aspects in engineering, technology and science concerning friction have impact on our 
daily lives [1|. As such frictional motion has been studied for centuries, but a complete physical 
understanding of friction is still lacking. For instance, the transition from stick to slip (the onset 
of dynamic sliding) is not well understood. Nevertheless, the initiation of dynamic sliding is an 
important aspect in many areas of science including fracture mechanics |]2||3| and seismology ||4]-[6|. 



The onset of dynamic sliding is often globally perceived as a uniform transition from sticking to 
sliding. In reality, however, it is a much more complex phenomenon. The shear stress distribution 
at an interface is generally nonuniform and reaches therefore the shear strength only at a narrow 
zone from which it might cause interface rupture. The repetition of such local slip events results 
in global sliding and provides a possible explanation of stick-slip behavior that is consistent 
with recent experiments, which showed that global sliding is preceded by local slip propagating 
over parts of the contact interface ||7l|8l. As shown in ||8l, these repeating precursors increase 
continuously their zone of propagation until a last precursor breaks the entire interface and 
causes global sliding. The propagation speed of interface ruptures was observed to range from 
slow P\^ to supersonic [3]. Moreover, the front speed of a single slip event can change along 



the propagation path Uj. By studying the stress field close to the interface, Ben-David et al. |10| 
observed experimentally that the rupture velocity of the detachment front is coupled to the local 
ratio of shear stress Ts to normal stress Og measured before slip initiation. 

Recently, numerical investigations 1 11-151 reproduced the general features of the experimental 
results of p 10| using simple spring-block models. In this Letter we study numerically the 
initial slip event using a finite element (FE) method (see also [16]), allowing us to access detailed 
information on the onset of dynamical sliding and to re-examine the hypothesis of Ben-David et 



al. 1 10 1 on the correlation between the slip front speed and Ts/cTs. The advantage of the FE method 
over the above-mentioned discrete techniques is the ability to reproduce correctly the mechanical 
behavior of continua (e.g. isotropy, elasticity). 

Model set-up 

The two-dimensional system under consideration consists of a rectangular isotropic elastic plate 
{w = 200 mm, h = 100 mm) in contact with a rigid plane [see Fig.jTTa)]. The corners of the plate are 
rounded to avoid stress singularities at the edges. To study this system we use a FE method with 
an explicit Newmark-j6 integration scheme 1 17| in plane stress incorporating an energy conserving 
contact algorithm. The material properties [see Table [T| correspond to polymethylmethacrylate 



glass (PMMA) which was also used in the experiments [ilOJ. We employ Rayleigh damping [18 19 1 
with mass and stiffness proportionality coefficient of and 0.1 jj.s, respectively. The deformable 
solid is discretized by regular quadrilateral elements (with element side ranging for different 
meshes from 0.67 mm to 2 mm) interpolating the displacement field linearly. 

A linearly distributed vertical displacement (u^ = 0.37 mm, u^ = 0.037 mm) is imposed at the top 
of the plate [see Fig.jTta)]. This loading is, after reaching equilibrium, complemented by applying a 
uniform horizontal velocity Vx = 10"^Cl, where Cl is the longitudinal wave speed in the deformable 
solid. The small value of the applied velocity ensures quasi-static loading conditions, similar to 
the experiments 1 10 1. The resulting stress distribution at the interface is nonuniform. Fig.jTTb) is a 



schematic depiction of the ratio of the local tangential traction tg to the contact pressure pg. These 
tractions (denoted with a subscript s) are measured at the moment preceding interface rupture 
and are referred to hereafter as static. The imposed loading conditions ensure a spontaneous 
nucleation of the first slip event inside the contact interface far from the edges [circle in Fig. [T|a)], 
because this is where the non-symmetric stress distribution reaches a critical value tg/ps > Ps, see 



Fig.ITTb). In the stick state, the tangential resistance of the interface is assumed to be proportional 
to the contact pressure p with a coefficient /.is- As for the slip state, this coefficient of proportionality 
[I is determined by the velocity (v) weakening friction law [see Fig. ITTc)] 

^ = Hk + il^s- l^k) exp (-|z;| V(i"s - l^k)/oc) , (1) 

which ensures a smooth transition from the static /,is to the kinetic /i^ friction coefficient governed by 
the transition parameter a. The parameters of the friction law are as well presented in Table [l| The 



local /is corresponds to experimental results as reported in Fig. 4(a-b) in 1 20 1 and is considerably 



higher than the global static coefficient of friction. An effect that was also observed in spring- 



block simulations 1 12 21 1. The local kinetic friction coefficient as well as the transition parameter 
were not measured in the experiments. Therefore, they were studied here qualitatively (see first 
paragraph of the following section) and eventually chosen arbitrarily. When the ratio of the local 
tangential traction to the contact pressure exceeds the static friction threshold {ts/ps > /is)/ slip 
occurs and propagates in one or both directions along the frictional interface. The dynamics of 
the slip fronts are determined by the parameters of the friction law (Eq. [l]) as well as by the local 
stress state. 

Results & Discussion 

We have conducted several simulations (not all presented in this Letter) and have observed 
different types of slip: crack-like (the entire interface between the crack tips is slipping), pulse-like 
(the slip region propagates along the interface within a narrow pulse) and mixed modes when a 
crack converts to pulses and vice versa. The propagation speed of the slip tip Vtip is related to the 
local stress state and seems not to depend on the type of slip. By studying the influence of the 
friction law parameters, we have observed that for an increasing (decreasing) difference between 
the static and the kinetic friction coefficients A/i = /is - /i)t/ the slip type tends to be crack-like 



(pulse-like). A higher transition parameter a causes slower slip propagation especially during 
slip initiation and slip arrest. 

In order to compare our numerical results with the experimental observations of Ben-David et 



al. 1 10 1, we present Kip as a function of the ratio of shear to normal stress measured before slip 
initiation. Here, the slip tip speed V,i^ is normalized to the longitudinal wave speed Cl and the 
local stress ratio is replaced by the local static ratio of tangential surface traction ts to contact 



pressure ps [see Fig.|2|. Our results confirm the experimentally 1 10 1 and numerically 1 14 1 observed 



general trend that the rupture propagation is faster for higher tg/ps ratios. For friction parameters 
/is = 1-3, jUfc = 1-0 and a = 0.1 m^/s^ the slip front velocities are in good quantitative agreement 



with the experimental results. Consistently with experiments 1 10 \, for the given type of loading 
(only at the top face), we do not observe slow fronts. Interestingly, we note that the rupture 
propagates considerably slower in the direction of the imposed shear load than in the opposite 
direction [compare solid with dashed line in Fig. |2l. These differences have not been reported in 
the experiments. 

To enable the separation of effects due to slip directionality and any other sources that might cause 
a non-unique relation between the tg/ps ratio and the rupture propagation speed we consider 
two additional simulations [Fig. |3l, where slip events are triggered at the edges. In order to 
increase the propagation distance [in comparison to Fig. |2| where the rupture propagating in the 
opposite direction of the imposed shear load arrests not far from the initiation zone] the kinetic 
coefficient of friction is reduced resulting in the following set of friction parameters: /Jg = 1-3/ 
[j.k = 0.6 and a = 0.1 m^s^. In all three cases the loading history of the body is identical up to the 
moment the tangential surface traction reaches the friction threshold, i.e., the initial stress state 
is the same for all simulations [see solid line in Fig. ffl. The slip propagation is then triggered 
by manually increasing the local tangential surface traction within small nucleation zones at the 



edges [Fig. |3ta),(c)]. Otherwise if the global shear load is slightly increased, rupture nucleates 
spontaneously far from the edges as before [Fig.lSpb)]. In case of spontaneous initiation [Fig.|3][d), 
solid line] the rupture propagates fast toward the edges and its velocity decreases along the path 
with a decreasing ratio tg/ps- Note that under some conditions we observe supersonic slip fronts, 
which were not observed in [10]. However, our results are consistent with rupture in bi-material 
interfaces where the stiffer material limits the propagation speed as observed experimentally and 



numerically by Coker et al. 1 22 1. For the two edge-triggered ruptures the slip propagates relatively 
slowly in the first phase, accelerates, reaches a maximum value (for maximal ratio tg/ps) arid 
decelerates afterwards [see Fig. |3|^d), dashed and dashed-dotted lines]. Although the triggered 
ruptures are unidirectional, there is no unique slip tip speed associated with a given fg/ps value. 
The maximal rupture velocity of the left-triggered slip does not exceed 60% of the maximal speed 
for the other two cases. 

As seen most clearly in Fig. ISla) the slip front (marked by a small white triangle) propagating 
at super-shear velocity follows the longitudinal wave (the circular white zone furthest from the 
nucleation zone), which modifies the local stress state at the interface. Therefore, looking at the 
dynamic ratio t^lpd measured in front of the slip tip, instead of examining the static ratio fg/Ps/ 
would allow to account for the dynamic nature of the slip propagation. 

Here, the location of the slip tip is determined to coincide with the position of the sticking node 
in front of the slipping nodes [see inset in Fig.H. According to this definition, the position of the 
rupture tip changes abruptly when the front advances. However, its velocity is computed in a 
continuous way as V^tip = I* I M, where I* is a characteristic distance (here I* = 0.67 mm) and A^ is 
the time interval that the rupture needs to advance this distance. 

In the context of discrete contact, we propose to analyze an instantaneous dynamic stress state 



(frf and pd) at the slip tip right after it jumps to a new position [see inset in Fig. |4|. The dynamic 
ratio td/pd differs significantly from the static one [compared in Fig. H, being changed by the 
longitudinal wave often preceding the slip front. It is worth noting that the value of the dynamic 
ratio is far from the critical value jj.s for a large part of the propagation path, which implies the 
need for a strong change of the stress state at the rupture tip within a short time. 

The relation between the velocity of the slip front and the dynamic ratio td/pd is depicted in Fig.|5j 
Compared to Fig. |3| the rupture triggered on the left is in better agreement with the other two 
(faster) slip fronts. Particularly, the slopes are more consistent for all curves and the range of 
velocities is smaller for a given ratio td/pd- Again it is confirmed that the character of the slip 
propagation is directionality dependent. For a given ratio td/pd, the slip fronts propagating in the 
direction opposite to the sliding are faster than the oncoming fronts [in Fig. |5| e.g., compare the 
dashed with the dashed-dotted curves]. Nonetheless, the difference between the curves cannot 
be only attributed to the directionality [in Fig. Ei note the two branches of the dashed and dashed- 
dotted curves]. The accelerating slip fronts show a faster rupture velocity than the decelerating 
ones for the same given ration tp/pd. Further, the general trend of faster rupture for higher t/p 
is lost [enclosed by the large circle in Fig. |5l; at a certain moment, the rupture speed starts to 
decrease rapidly with increasing td/pd along the propagation path. We observe this phenomenon 
only for slip fronts advancing against the sliding direction. Regardless of the simplicity of the 
static criterion ts/ps and the consistency of the dynamic criterion td/pd, a stress ratio does not seem 
able to provide a fully reliable estimation of the velocity of the slip propagation. 

The lack of generality of the velocity criteria based on the ratio of the tangential traction to the 
contact pressure t/p suggests an independent consideration of t and p. It was proposed pO| that 
the propagation of the slip front is related to the energy densities Ug, stored at the front tip, and Ur, 



needed to advance the slip front. We propose a heuristic energy density at the contact interface as 

U(p,t) = (2{l + v)f + p^)/2E. (2) 

The density of stored energy Us = U (pd, t^) is measured locally at the slip tip at the moment 
the front advances one length parameter I*, similarly to the dynamic ratio td/pd- The density 
of rupture energy U,- = U (p,-, fJ-sPr) is computed at the same material point just before the front 
advances another l*, i.e., when the ratio of tangential traction to contact pressure reaches the static 
coefficient of friction {tr/p,- = /^s) [see inset in Fig.[6l. 



The normalized rupture velocity is depicted in Fig. |6] as a function of the change of the energy 
density at the slip tip ALT = U,- - Us normalized by the stored energy density Us. The data of all 
three cases collapse within a narrow region properly described by 

VJc^ = a + b exp(-c ^JAU/Us) , (3) 



where a, b and c are fitting parameters [see Fig. [6|. No differences due to the directionality of the 
slip propagation nor any other reason that caused branching for the previously studied criteria 
are now present. This shows that the energy density criterion is able to account for the dynamics 
of slip events at bi-material interfaces. Note that tails of data points falling outside of the fit range 
occur when the slip fronts start to decelerate rapidly before arresting. 

Conclusion 

In this Letter it is demonstrated that the static ratio of shear to normal stress |T0|[T4| is not 
a sufficient criterion for determining the speed of slip fronts. The use of the dynamic ratio, 
measured in front of the slip tip, improves the estimation of this speed. However, for our set- 
up we observed that, given a stress ratio (static or dynamic), the front going in the direction 
of the sliding is always slower than the front propagating in the opposite direction. Moreover, 

8 



the decelerating fronts are also slower than the accelerating ones. The energetic criterion we 
propose eliminates these effects and highlights the similarities between the rupture of frictional 



interfaces [3] and crack propagation |22|. It is hoped that these findings motivate experimental 
work to access dynamic stress field measurements as well as theoretical studies to extend the 
principles of fracture mechanics to problems of frictional sliding. 
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Table 1 Friction and material parameters corresponding to polymethylmethacrylate glass, PMMA 



Parameter 



Material: 
Young's modulus E 
Poisson's ratio v 
Density p 

Longitudinal wave speed Cl 
Transverse wave speed Cs 

Friction: 
Static friction coefficient /is 
Kinetic friction coefficient [J.k 
Transition parameter a 



2.6 GPa 

0.37 

1200 kg/m^ 

1584 m/s 

890 m/s 

1.3 
0.6; 1.0 

0.1 mVs^ 




y /J // / ^ Rt7 



\i, 



CD 
CJ 




h % 

I o 



CJ 








K ^s = 1.3 (c) 
\ \ A^it = 0.6 

V 


« = 0.005 

tt = 0.100 

a = 0.500 


mys^ 

mys^ 



12 3 

slip velocity \v\ (m/s) 



Fig. 1 Two-dimensional set-up of the problem: (a) a thin rectangular plate in contact with a rigid plane 
is loaded on the top by a linearly distributed imposed displacement Uy and a uniform velocity Vx) (b) the 
nonuniform distribution of shear to normal tractions (fg and ps/ respectively) at the interface causes a first 
slip nucleation far from the edges [fs/ps > /^s marked by a circle in (a)]; (c) the change of the friction coeffi- 
cient with respect to the material slip velocity v is governed by the parameter a [see Eq. Tl 
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Fig. 2 Comparison of numerical results with experimental observations by Ben-David et al. [10|. The nor- 
malized rupture velocity is reported with respect to the static ratio of local tangential surface traction ts to 
contact pressure ps- Friction parameters are /Js = 1-3, /jjt = 1-0 and a - 0.1 m^s^ 
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Fig. 3 Three different sHp events are presented for the same initial stress state (before triggering or spon- 
taneous initiation), h:\stantaneous material velocity is shown for the slip event (a) triggered at the left 
edge, (b) spontaneously initiated far from the edges and (c) triggered at the right edge. Colors from blue 
to white denote material velocities ranging from Om/s to 2m/s, respectively. The starting point of each 
event is marked with a square, a circle and a triangle, respectively. Small white triangles show the lo- 
cation of the tip of the slip front. Black arrows indicate the direction of the imposed global shear load, 
whereas white arrows show the direction of the rupture propagation, (d) The normalized rupture ve- 
locity for all three cases is depicted with respect to the local static ratio of tangential traction tg to contact 
pressure p, (data close to the triggering zone are not shown). Friction parameters are /.is = 1-3, ^k - 0.6 
and a = 0.1 m^/s^ 



13 








50 



100 

X (mm) 



150 



200 



Fig. 4 For each point along the interface x the dynamic stress ratio f^f/p^ is plotted at the moment when 
the slip front arrives at this location. Note that contrary to the reported static stress ratio tgjps, this is not 
an instantaneous picture but an assembly of results over the entire time of propagation. Data close to the 
triggering zones are omitted. (Inset) The dynamic values t^^ and pd are measured at the sticking node in 
front of the slipping region at the moment f"*"^ when the previous node starts to slip 
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Fig. 5 The normalized rupture velocity is plotted with respect to the dynamic traction ratio frf/p^^ for all 
three slip events. Data close to the triggering zones are omitted 
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Fig. 6 The normahzed rupture velocity is depicted with respect to a dynamic criterion based on a heuris- 
tic surface energy density M. The ratio (LT,- - LTs) /iis represents the proportion of the energy change 
Mr - Us at the slip tip needed to advance the rupture front with respect to the locally stored energy Us- 
The gray area is a data fit based on Eq.|3]with a - 0.76 + 0.07, b - 1.80 and c - 3.05. (Inset) Us and U,- are 
measured when the slip tip, respectively, reaches the observation point and overpasses it 
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